Condensation Phenomena in Nano-Pores - a Monte Carlo Study 



Raja Paul 1 , Heiko Rieger 1 
1 Theoretische Physik, Universitdt des Saarlandes, 66041 Saarbrilcken, Germany. 

The non-equilibrium dynamics of condensation phenomena in nano-pores is studied via Monte 
Carlo simulation of a lattice gas model. Hysteretic behavior of the particle density as a function of 
the density of a reservoir is obtained for various pore geometries in two and three dimensions. The 
shape of the hysteresis loops depend on the characteristics of the pore geometry. The evaporation 
of particles from a pore can be fitted to a stretched exponential decay of the particle density pf(t) ~ 
exp [— (t/r)^] . Phase separation dynamics inside the pore is effectively described by a random walk 
of the non- wetting phases. Domain evolution is significantly slowed down in presence of random 
wall-particle potential and gives rise to a temperature dependent growth exponent. On the other 
hand roughness of the pore wall only delays the onset of a pure domain growth. 

PACS numbers: 64.60-i, 81.07-b, 47.55. Mh, 75.40. Gb 



I. INTRODUCTION 



Adsorption and desorption isotherms of a gas con- 
densed into nano pores show hysteresis and have become 
very a useful tool for the classification of porous mate- 
rials p], 0, 0. Due to the effect of surface tension the 
formation of meniscii inside the pore is delayed on the 
desorption branch, resulting in a non-vanishing hystere- 
sis area of the sorption isotherms. A theoretical analysis 
and detailed description of the early research on this phe- 
nomenon is described in a review by Everett Q. Recent 
investigation from a point of view of the stability of ad- 
sorbed multilayers 5] and other 0, 0, Q analytical and 
numerical approaches together with density functional 
theory also proposed that the hysteresis phenomenon is 
an intrinsic property of the phase transition in a single 
idealized pore and arises from the existence of metastable 
states. 

Hysteresis in real materials is the collective phenomena 
involving interconnected network of simple pores 0, [HI • 
However , qu ite recently it has been shown by numerical 
analysis fTH and experiments 0, 0, ^( that hysteresis 
can occur in isolated pores also. Moreover the shape of 
the hysteresis loops are influenced by the characteristic 
features of the pore geometry. In this paper, using exten- 
sive Monte Carlo simulations, we intend to characterize 
the sorption isotherms for different pore geometries and 
compare our results directly to the experimental obser- 
vations. 

Furthermore, we study the phase separation kinetics 
of a binary liquid in nano-pores at low temperatures. 
Previous work [15|, jig, ll7| showed that a binary liquid, 
unlike Ising-like complete phase separation, does not sep- 
arate into two phases completely. Instead the adsorbed 
material forms many small domains far below the co- 
existence region. An explanation for this behavior has 
been suggested on the basis of random-field Ising model 
pH [HI [23 • On the contrary we found in case of a single- 
pore model with no randomness, confinement in a small 
pore slows down domain growth in certain regions of the 
wetting phase diagram [2l| and as a result macroscopic 



phase separation is not observed on short time scales. 
Moreover the late stage domain evolution, obtained from 
a two point correlation function C(r, t) follows the same 
Lifshitz-Slvozov [23 growth law of t 1 / 3 . Unlike earlier 
research [23| of domain evolution in porous networks, we 
concentrate on a single pore which is more applicable for 
systems with low porosity (Si mesopores [l^). Moreover, 
randomness in terms of irregular pore structure or pres- 
ence of impurity atoms in the pore wall, which is inherent 
in real systems, may have drastic effect on the domain 
growth. 

In this paper we focus on a single pore. An average 
over many isolated pores corresponding to an assembly 
of non-interconnected (unlike Vycor) pores can be exper- 
imentally developed [l4| out of a B-doped Si wafer via 
electrochemical etching. The structure of the rest of the 
paper is as follows. First we overview, in section II, the 
theoretical model and the Monte Carlo technique that 
we use in our simulation. In section III A, we discuss the 
Hysteresis phenomenon for different pore geometries in 2 
and 3 dimensions and discuss how one can differentiate 
the pore structure from the hysteresis loops. In section 
III B, we focus on the evaporation of particles from the 
pore and discuss how the density in the pore reaches a 
new equilibrium value. Sections III C, D are entirely de- 
voted to the domain evolution in the pore environment 
where we propose a random-walk picture of the domain 
growth above a critical size comparable to the pore di- 
ameter. Then we perform a quantitative study of growth 
phenomenon for a simple pore as well as for pores with 
defects. Finally section IV finishes the paper with a dis- 
cussion. 



II. SIMULATION MODEL 



A standard model for a binary liquid mixture, is the 
Ising lattice-gas model with spin occupancy variables &i 
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= 0, 1 governed by the Hamiltonian: 
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where W pp and W^p are the particle-particle and wall- 
particle couplings respectively. We denote an occupied 
site as a particle and an empty site as a vacancy. Exper- 
imental observations suggest that in case of a glass(Si) 
adsorbent, the pore wall has a very strong affinity to- 
wards the adsorbed gases, which implies W wp > W pp . 
Most of our simulations, unless otherwise specified, are 
performed at a fixed "wettability" W wp /W pp = 1.5. A 
change in this value does not qualitatively affect our re- 
sults as long as W wp > W pp . The geometry of the pores 
is chosen to be a rectangle of size L x h in two dimension 
and parallelepiped of size L x h x h in three dimension 
with h <C L. Standard conserved order parameter dy- 
namics (Kawasaki) is employed to study the diffusion and 
the domain growth kinetics in the nano pores. Nearest 
neighbor spins and crj are exchanged with Metropolis 
acceptance probability 



P({a} — > {a'}) = min 



1 exp{-PHW}) 
X ' exp(-(3H{a}) 



(2) 



where a and a' represents the old and new spin configu- 
ration respectively. To improve the speed of simulations 
we always choose a particle with unit probability and de- 
fine 1 MC step containing a number of trial-updates that 
is equal to the number of particles present in the system. 



III. NUMERICAL RESULTS 

A. Hysteresis 

To study hysteresis in the adsorption/desorption pro- 
cess within a pore we attach one or two reservoirs depend- 
ing upon the geometry of the pore. To avoid asymmetric 
diffusion into the pore, no periodic boundary condition 
is applied between the reservoirs. The particle density in 
the reservoirs is kept constant by adding (removing) parti- 
cle at a randomly chosen position in the reservoir as soon 
as a particle diffuses into the pore (reservoir). A snapshot 
of the hysteresis phenomenon for a simple pore geometry 
in 2d is shown in Fig. Initially, both the pore and the 
reservoirs are kept empty. The density of particles in the 
reservoirs(p res ) are then slowly increased. Particles are 
immediately adsorbed and form a single layer along the 
pore walls. At this stage the density in the pores(p/), 
rises sharply to a non-zero value which corresponds to 
the first jump of adsorption isotherms (Fig. |5J). Then pf 
forms a long plateau until two semicircular (or hemispher- 
ical in 3d) meniscii(a meniscus for one-end open pore) are 
formed somewhere in the middle of the pore at a high 
reservoir density and move apart from each other to fill 
the entire pore; a second jump in pf is observed. We 
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FIG. 1: Snapshots of the hysteresis in a simple nano-pore 
of size L — 100 and h = 11 at temperature T — 0.3. The 
shaded regions and the black spots on both sides of the pore 
are respectively the attached reservoirs and the particles in 
it. For each value of p res the system is equilibrated up to 
t = 2 22 . Snapshots (a) correspond to the adsorption and (b) 
correspond to the desorption isotherms. 



call this complete pore filling and further increase in p res 
does not change pf. 

Next, we slowly decrease the reservoir density and 
record the corresponding filling fraction pf after the sys- 
tem is equilibrated. On desorption, however at a much 
smaller reservoir density p res ^ we see the semicircular 
meniscii to form on the open ends of the pore, followed 
by a sharp fall in the pore density p/, as the meniscii 
approaches towards each other. Finally the desorption 
curve follows the same path as that of the adsorption. 
Since the diffusion occurs very slowly into the pores, a 
long waiting time t w is needed to equilibrate the whole 
system for each value of p re s- 

Our simulation results are shown in Fig. where col- 
umn (a) and (b) shows the hysteresis in 2 and 3 dimen- 
sions respectively. The simulations are performed with 
systems of length L = 100 and 50 at temperatures T = 
0.3, 0.6 in 2d and 3d respectively. Reservoir lengths are 
kept fixed at L r = 10. Ink-bottle pores are equally di- 
vided into two parts; one with the narrow tube of diame- 
ter hi and the other is relatively wider sack with diameter 
/i2- To make sure that pore density pf does not change 
any more, each data point of the sorption-isotherms is 
equilibrated for t w = 2 22 time steps. Finally for each 
hysteresis loop data are averaged over 50 ensembles. It 
is evident from the sorption branches of II-IIIa,b, that 
they are much steeper than in Ia,b. This is due the for- 
mation of two meniscii on both sides of the pore, which 
accelerate the filling and emptying procedure. Such char- 
acteristic feature distinguishes between a one-end open 
and both-end open pore. Further, an increase in pore di- 
ameter requires a higher reservoir density to initiate the 
formation of meniscii, which effectively delays the onset 
of adsorption saturation. As a result the the loop area 
for bigger pore is also increased as shown in IIa,b and 
IIIa,b. This features can be used to compare pores of 
different sizes. 
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FIG. 2: Adsorption (filled circles) and desorpt ion (empty cir- 
cles) isotherms for nano pores, (a): Filling fraction pf of the 
pore as a function of the reservoir density p res in 2 dimension. 
The equilibration time t w for each data point 2 22 . Ia has one 
reservoir of length 10 attached to the open end. Ha, Ilia and 
IVa has 2 reservoir attached on both sides, (b): Hysteresis 
isotherms in 3 dimension. Equilibration time tw — z , pore 
length L — 50 and reservoir length is 10. For each hysteresis 
isotherm data are averaged over 50 disorder samples. 



In experiments one observes a two step decrease 
of the desorption isotherms corresponding to ink-bottle 
pores. Our data, as shown in figure Vb, for a both-end 
open pore in 3d of similar geometry, agrees quite well with 
the experimental predictions. The two step decrease of 
the desorption isotherm arises from the sequential emp- 
tying of the different sections (of diameter hi and Z^) of 
the pore. The heights of the two steps suggest the wider 
section of the pore emptied earlier than the smaller sec- 
tion. However, ink-bottle pores in figure IVa,b(2<i-both 
end open) and Va(3<i-one end open) do not exhibit such 
special characteristics. Absence of this feature may be 



due to a small temperature or a small system size with 
insufficient ^2/^1 ratio. 



B. Evaporation 

If a partially or completely filled pore is kept in vac- 
uum, the density inside it decreases with time. This is 
what we call evaporation and investigate it for different 
pore geometries. In experiments 14] one measures the 
vapor pressure change inside a previously equilibrated 
pore subject to the pressure variation in the reservoir. 
However, we carry out the simulation in a slightly dif- 
ferent way by keeping the initial pore density just above 
the desorption threshold. 

To study the variation in pore density pf as a function 
of time we fill pores completely and allow it to evapo- 
rate in an empty reservoir (vacuum). The change in pf 
has been recorded as a function of time t. Decay of the 
pore density can be very well described by a stretched 
exponential law, 



p f (t) ~ exp [-(t/rf 



(3) 



The simulation is carried out in pores with both sim- 
ple and ink-bottle geometry. For simple pore we study 
the evaporation at temperature T =0.4, for systems with 
L = 128, 256, 400, 512 and h =7 and finally average 
over 50 ensembles. The pore density pf as a function of 
time t is plotted in Fig. 01(a) in a log-linear scale. It is 
noticed that a pure stretched exponential decay is found 
only above a certain value of the pore filling fraction as 
shown by continuous lines. These values of pf are noth- 
ing but the ratio between the number of surface to bulk 
molecules and decreases as L becomes larger. The surface 
molecules, are attached rather strongly (W wp > W pp ) to 
the walls and slow down the evaporation rate for a while 
but finally drop off suddenly to zero. As the length of 
the system is increased, the particles deep inside the pore 
requires much more trial attempts to diffuse till the open 
end which effectively increases the evaporation time for 
longer pores. This leads to the decay exponent f3 to de- 
crease with L, as shown in the inset of (a). 

A similar study in case of ink-bottle pores was carried 
out with systems of fixed L = 200 and hi = 7 for differ- 
ent values of /12 = 7, 21, 35, 49 as shown in Fig. 01(b). 
The temperature is kept fixed at T = 0.4 and data are 
averaged over 50 ensembles. Like the previous case, here 
also, pf shows a pure stretched exponential decay above 
a certain value, which becomes smaller as h^ is increased. 
Moreover the effective length(L + /i2) of the pore increase 
with h2 results in the decay exponent (3 (inset of (b)) to 
drop off similarly to the simple pore. 

A further analysis (not shown) on both simple and ink- 
bottle pore shows that f3 is independent of the tempera- 
ture. 
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FIG. 3: (a): Evaporation from simple pores of diameter h = 7 
and different lengths L =128, 256, 400 and 512 from left to 
right at temperature T = 0.4. Data for the filling fraction 
Pf(t) are shown in a log- linear scale - note that the decay 
is initially slower than exponential. The full line connects 
the data points that can very well be fitted to a stretched 
exponential with the exponent (3 given in the inset. Inset: 
The decay exponent /3 as a function of the pore length L. 
(b): pf(t) for ink-bottle pores of fixed length L — 200, tube 
diameter hi = 7 at T = 0.4 for different sack width hi =7, 
21, 25, 49 from left to right. The decay law is similar to the 
simple pores. Inset: (3 as a function of hi. 




C. Domain evolution and Random Walks 

In this section we study the temporal evolution of do- 
main structures inside pores starting with a random(high 
temperature) initial configuration of given density. Now 
we use periodic boundary conditions between the two 
ends of the pore and the total number of particle is con- 
served. Fig. Ufa), shows the domain structure-profile 
at different time steps when the system is quenched to 
T = 0.3 from T = oo. Fig. Q|b) is the snap-shot of a 
cut along the pore axis. The black regions correspond 
to the occupied sites and are called "particle-domains", 
whereas empty white regions are termed "blobs" . 



FIG. 4: (a): Domain evolution in nano-pores. Figures are 
produced by taking snap-shot for a system of length L = 
300 and height h = 7 at temperature T = 0.3. (b): Axial 
snapshot of the domains for the same system in Fig. (a). A 
horizontal cut through the figure at any instant t will give the 
location and size of domains/blobs along the axis of the pore. 
The time is plotted in logarithmic scale, (c): Time evolution 
of a single rectangular blob of linear size lb = 13 in a pore of 
dimension 200 x 7 at T = 0.3. 
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In the initial stage, the growth is dominated by nu- 
cleation and spinodal decomposition. As soon as the 
domain- or the blob-size becomes comparable to the di- 
ameter of the pore, the above two mechanisms do not 
work any more. Transverse directional growth is com- 
pletely stopped because of the pore wall and the horizon- 
tal movement of the blobs is slowed down by the presence 
of the particle-domains. 

The snapshots in Fig. 0] confirm that the random mo- 
tion of the blobs plays an important role in the late 
stage of growth. A closer look to Fig. EJb), which is 
horizontally the occupancy along the pore-axis and ver- 
tically time in log scale, shows in the late stage, the 
blobs move to and fro along the axis of the pore and 
penetrate through the particle-domains to coalesce with 
the neighbors. During this process it also transfers a 
holes (vacancies, white regions) from its surface to the 
neighboring blobs. 

To elucidate the random motion, we study the time 
evolution of a single blob as shown in FigOJc). Initially 
it is a perfect rectangle of linear size(Z&) comparable to 
the pore diameter (ft) and placed in the center of the pore. 
The temperature is kept fixed at T = 0.3. The system is 
allowed to evolve and the mean square deviation (x 2 (t)) 
of the center of mass (CAM) of the blob along the axis of 
the pore is measured at each time step. Since the bound- 
ary of the blob perpendicular to the wall, fluctuates very 
rapidly, the true COM may not lie on the geometrical 
axis of the pore. In this case we trace the actual COM 
and project it on to the axis. The size of the blob has 
to be chosen large enough to avoid disintegration of the 
main blob. 

Since we start with an exactly rectangular blob, a true 
random- walk motion is not immediately observed. In the 
early stages, the blob walls, perpendicular to the pore 
axis roughens, leading to (x 2 (t)) ~ t 2//3 as in surface 
roughening described by the KPZ-equation |24|. At late 
stages (t > 10 4 ), the blob performs a random walk with 
a blob-size dependent diffusion constant £>(/&), 

(x 2 (t))=V(l b )t. (4) 

Fig-Efa), shows the mean square displacement (x 2 (t)) 
of the COM of blobs with l b =11, 15, 19, 23, 27, 31, 
35, 39 as a function of time t. For each values of £& the 
simulation is carried out in a system with fixed L = 200, 
ft = 7 at temperature T = 0.3 and finally averaged over 
1000 ensembles. It is evident from the figure that the 
onset of a true random- walk is delayed for bigger blobs. 
A scaling form for both small and large time scale regime 
is, 

(x 2 (l b ,t))<XT 6 (l b ) f(t/r(l b ) (5) 

with 

-/ x f x e for x <C 1 fn x 
x for x»l. (6) 
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FIG. 5: (a): Mean-square deviation of blobs of size h = 11, 
15, 19, 23, 27, 31, 35, 39 for a system size L = 200, ft = 7 and 
the temperature T = 0.3. (b): Shows the power low decay 
of the diffusion constant T>(h) for different blob sizes. The 
exponent is found to be 7 = 0.68. 

Assuming, r{l b ) ~ h z •> from Eq.0 one can write for large 
time t ^> r, 

(x 2 (l b ,t))<xl b z{e -V t, (7) 
which readily gives, 

v{i b ) ~ i; z{1 - 9) = h' 1 . (8) 

The scaling form in Eq. [51 gives quite reasonable data 
collapse for = 2/3 and z = 2, as shown in Fig. \^a). 
Substituting these value in Eq. [S| we obtain the value of 
7 = 2/3. Diffusion constant V evaluated from the asymp- 
totic behavior [c. f. Fig[Sfa)] is plotted in FigOb). From 
the slope of V vs. curve in log-log scale we estimate 
the 7 = 0.68 which agrees quite well with the previous 
value 2/3, obtained from the scaling. Note that the re- 
lation V oc 1 / deviates from the naive expectation 
V ex ksTrj/lb [2l|, |25j, where T, r] are the Boltz- 
mann' constant, temperature and viscosity coefficient of 
the medium respectively. 
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D. Correlation function and Domain growth 
exponent 

An alternative way to study domain growth is via the 
measurement of the correlation function C(r, t) along the 
axis of the nano-pore in a similar fashion as in the Ising 
model. 



(a) 



C(r,t) = (S(0,t)S(r,t)) - (Sf 



(9) 



where S(r, t) = 2<j(r, t) — 1, the Lattice gas- variable, takes 
the values -1, 1 for a = 0, 1 respectively. Due to this 
transformation, C(r, t) falls off with r in an oscillatory 
fashion, as is shown in the inset of Fig. EJa). Following 
Huse 27], we also define the length scale or the typical 
domain size R(t) of the system as the position of the first 
zero of C(r, t). Calculating C(r, t) using Eq.Elwe extract 
R(t) by fitting the three or four points in C(r, t) closest 
to its first zero to a quadratic function of r and defining 
R(t) as the value of r for which the function vanishes. 
At very late stage (t > r(T)), the asymptotic domain size 
R(t) grows as a power-law, 



R(t) ~ t 



(10) 



where 1/iJj is the growth exponent. The onset of such a 
power law regime shifts towards the early time scales as 
the temperature is increased. 

We studied the domain growth for both (i) simple and 
(ii) complex pores. By "simple" we specify a pore with 
no geometrical defects in the wall and has a fixed value of 
W wp . "Complex" pores are further subcatagorized into 
(ii)-a: simple pore with random W wp G [0, 1.5] and (ii)- 
b: pore with fixed W wp but geometrical defects along the 
wall. Let us discuss about the "simple" pores first. The 
domain evolution in this case studied with a system of 
length L = 500 and height h = 7 at different temperature 
T = 0.25, 0.3, 0.35, 0.4. For each temperature we aver- 
aged over 300 ensembles. The average length scale or the 
domain size R(t) as a function of the time t is plotted in 
the Fig.UJa). The growth exponents 1/^, extracted from 
this figure, for different values of temperatures, is shown 
by upper curve in Fig.|^b). This exponent appears to be 
~ 1/3 independent of the temperature. On the other 
hand, for the case (ii)-a with defects in terms of random 
wall-particle coupling W wp G [0,1.5], we carried out a 
similar kind of study as described above. It is observed 
that the growth process is drastically slowed down by the 
pinning -effect 281 of non-wetting sites located randomly 
in the pore walls. And as a consequence the exponents 
are also reduced quite significantly and found to have a 
temperature dependence as well, as shown in the bottom 
of Fig.Etft). 

It is also noticed, for the case (ii)-b, that a periodic 
structure-defect of the pore wall has a strong influence 
on the domain growth. In Fig. Ufa) we demonstrate the 
effect of a variable periodicity Rl for system of length 
L = 512 with two fixed pore diameters h\ = 9, h<i = 15. 
The simulations are done at a constant temperature 
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FIG. 6: (a): Average domain size R(t) as a function of time 
t for different quenches to T = 0.4, 0.35, 0.3. The system size 
L = 500, h — 7 and the density pf = 0.7. Inset : Correlation 
function C(r, t) at various times t after a quench to T — 0.35. 
The domain size R{t) is defined as the distance of the first zero 
from the origin, (b): Growth exponent plotted against 
the temperature T. The upper curve, for fixed W wp , is almost 
a constant with T, whereas the lower one, for random W wp , 
shows a steep decrease for smaller values of T. 



T = 0.35 and for each value of Rl, the data are averaged 
over 300 ensembles. Domains grow quite fast in either 
part of the hi or diameters, but finally due to the 
bottle-neck effect it becomes extremely hard to transfer 
the particle-domains or blobs from one side of the narrow 
part (diameter hi) of the pore to the other. Consequently 
the growth process is slowed down. But in the late stage, 
when the domain size R(t) becomes comparable to the 
periodicity Rl, the blobs are in similar environment like 
in a simple pore and a true power law growth sets on. 
As a result a longer periodicity corresponds to a late 
initiation of the power law growth. We carried out the 
simulation for 4 different values of periodicity Rl =4, 8, 
16, 32 and for the first 3 values of Rl we found a late 
stage growth exponent = 0.33, but for Rl = 32, 
the power law growth in not yet initiated. For such a 
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large periodicity one certainly has to wait extremely 
long. Further, we carried out a similar study by varying 
the pore diameter hi while keeping the periodicity R L 
fixed. The numerical data are obtained with a system 
of L = 512, hi = 9, R L , T = 0.35 and h 2 =13, 15, 
17, 23. Finally for each value hi we average over 300 
ensembles. The domain size R(t) is plotted in Fig.^b). 
In the initial stage of the growth, as described above, 
the domain grows quite independently in either part of 
the pore with hi and hi diameters. A true power law 
growth, which starts after R(t) > Rl, is slightly delayed 
for larger hi. As described above, to attain a pore 
environment similar to the simple pore it is necessary to 
equilibrate the top and bottom of the wider part (hi) of 
the pore which takes longer as hi is increased. Finally, 
the growth exponent extracted in the late stage, found 
to have the same value = 0.33, as that of the pure 
model. 



IV. DISCUSSION 

Using extensive Monte Carlo simulations we have 
shown how hysteresis arises in nano-pores for differ- 
ent pore structures. The characteristics of the sorption 
branches are influenced by the shape and size of the pore 
geometry. Since the hysteresis in nano-pores occur due 
to diffusion, the temperature has to be chosen very care- 
fully to avoid slow diffusion rate at low temperatures or 
vanishing hysteresis- loop at high temperatures. Absence 
of two-step desorption branch in 2d ink-bottle pores may 
be due to the effect of a small temperate as found in 
the experiment Choice of pore diameters for one- 

end open ink-bottle geometry is also an important factor. 
Because of surface tension effect, a large curvature ratio 
7^ / 7^- > can cause a huge pressure difference between the 
narrow and wide parts of the pore. Increase in curvature 
corresponds to a decrease in vapor pressure; described by 
the Kelvin equation. As an effect if the narrow part of 
the pore is filled early, can block the particle to diffuse 
into the wider part, unless a very high density is reached 
in the reservoir. 

The stretched exponential decay of the pore density 
Pf(t) also agrees qualitatively with the experiments [l4| . 
One simple reason of such a behavior of pf(t) could be 
due to the formation of meniscus at the open end of the 
pore which eventually restricts the rate of evaporation 
depending upon the radius of its curvature. 

For the phase separation of binary liquids in the pore 
environment, model B j2<j corresponding to dynamics 
with conserved order parameter, is not well suited, as it 
does not account for the transport of the order parame- 
ter by hydrodynamic flow. Modifications to model B by 
adding an "advection" term describes the system quite 

weii [23. 

Late stage coarsening in pore system that we studied is 
effectively driven by two mechanisms: (l)Transfer of holes 
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FIG. 7: (a): Average domain size R(t) of a defected pore- 
wall for periodicities Rl — 4, 8, 16 and 32 at temperature 
T — 0.35. The system has fixed length L — 512 and diameters 
hi = 9 and hi — 15. Domain growth slows down for larger 
periodicity, (b): R(t) for different pore diameters hi =13, 15, 
17 and 23 at fixed L = 512, hi = 9, R L = 16, and T = 0.35. 
The domain growth remain constant while varying hi. 



from the surface of one blob to a neighboring blob and 
(2)the transfer of particles from one side of the blob to the 
other along the pore wall. Owing to the first mechanism a 
blob shrinks and the neighboring blobs grow, whereas due 
to the second one a blob moves to and fro as a whole and 
coalesce with another. Our numerical study for a single 
blob (section 111(C)), account for the contribution arises 
from (2) only and gives rise to the diffusion constant V ~ 

Ifr 2 ^ 3 of a blob of size The first mechanism modifies 
the single-blob diffusion constant in a non-trivial way, 
which is difficult to estimate. Finally, the superposition 
of random motion of blobs due to mechanism (2) and 
the hole transfer mechanism in (1) leads to a late-stage 
growth law R(t) ~ t 1 / 3 , independent of their individual 
contributions. 

The domain growth exponent is significantly de- 
creased in presence of random wall-particle potential 
W wp which in practice is due to the effect of impurity 
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atoms at the pore walls. The temperature dependence of 
the exponent may come from a logarithmic scaling [2^ 
of barrier energy of the domains. On the other hand geo- 
metrical defects appear to slow down the growth process 
only in the early time regime. But in the late stage, as 
the blob size becomes larger than the wave length of the 
geometrical disorder, which in this case is the periodic- 
ity of the wall-defect, the domain growth shows a pure 



behavior. 
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